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Abstract. Extreme objects such as X-ray binaries, AGN, or y-ray bursters harbor high 
energy plasmas whose properties are not well understood yet. Not only are they responsible 
of the hard X- and y-ray emission we observe but also they have a strong influence on the 
main dynamics and energetics of these objects themselves. 

Here we present a new kinetic code that solves the evolution equations for particles and pho- 
tons around compact objects. It produces spectra that will be compared with observations 
from Simbol-X to constrain the radiation and acceleration processes in these objects. 
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1. Introduction 

Investigating the time evolution of interact- 
ing particles and photons implies to deal with 
integro-differential equations that cannot be 
solved analytically, and numerical simulations 
are required. We have started the develop- 
ment of a new code that deals with parti- 
cles from sub- to ultra-relativistic regime, and 
includes interactions with the photon distri- 
bution such as self-absorbed synchrotron and 
Comptonization. It is time dependent and so 
will be able to address the variability observed 
in the sources light curves and spectra. 

2. The code 

The code solves simultaneously two identical 
equations that describe the time evolution of 
both the distribution of particles in the the mo- 



mentum space (x = p/mc) and that of photons 
in the energy space (x = hv/mc 2 ) : 

d,N = d x (A + d x (DN)) + a nj - 

1 esc 

The sink term N/T £SC represents the escape of 
particles/photons from the system. The source 
term ginj includes the direct injection of parti- 
cles/photons in the system (e.g. loading of mat- 
ter or seed photons from an accretion disc) as 
well as interactions between the two popula- 
tions (see hereafter). When possible/required, 
these interactions are treated in the Fokker- 
Planck approximation that leads to the first 
two terms corresponding to advection (A) (i.e. 
heating/cooling) and diffusion (D) in the mo- 
mentum/energy space. 

The distributions are discretized in bins. 
Since the Courant condition for explicit meth- 
ods sets a very small time step when the energy 
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range spans over several orders of magnitude, 
a semi-implicit method is used. The Chang- 



Cooper met hod (|Chang & Cooper|[1970| |Park 
|& Petr osian 19996) is inaccurate when solv- 
ing the FP equation with self-absorbed syn- 
chrotron for it leads to a poor energy conser- 
vation. Rather, we use a scheme based on the 
equations written above, which insures a num- 
ber and energy conservation to machine preci- 
sion. 

3. Radiation processes 

So far, two radiation processes have been im- 
plemented. 

Self-absorbed cyclo-synchrotron radiation: 

The power spectrum of a single electron is tab- 
ulated from a combination of asymptotic ex- 
pressions to get a good accuracy in all regimes 
(Ghise TIrri & Svensson|199l| ). Contributions of 
the cyclo-synchrotron radiation to the FP coef- 
ficients in the equation for the particles are cal- 
culated from expressions given in Ghis ellini et| 
|al~| ( fT998] l. 
Compton scattering: 

Including the Compton contribution for all en- 
ergy regimes requires to compute the exact 
distribution P(po, u>o; ai) of photons resulting 
from the interaction of particles of momentum 
Po with isotropic photons of energy wq- In the 
small angle scattering limit, the contribution 
to the Fokker-Planck coefficients is computed 
from the first moments of the scattered distri- 
bution. However, for the photon equation or for 
large angle scattering of particles, the FP ap- 
proximation is not relevant and an exact inte- 
gral treatment is used (see Nayaks hin & Melia] 
[T998] l. 




Fig. 1. Time evolution of the particle popula- 
tion when mono-energetic particles are con- 
stantly injected in an empty system (the time 
is normalized by the light crossing time of the 
system R/c and the distribution by Rcr T ). 

self-absorbed photons is much more efficient. 
Fig. |4| shows res ults very similar to those of 
Ghise llini et al. | ( |1998| : as the the corona is 
filled by particles, their distribution evolves to 
a Maxwellian one on a few synchrotron times 
scales, showing that exchange of photons is a 
very efficient thermalization process. 

5. Conclusion 

We have presented the basic properties of the 
code at its present stage of development and 
one first result. New features will be added in 
order to give a complete modelling of high en- 
ergy plasmas in microquasars, AGN, and y- 
ray bursts... Among them: Coulomb interac- 
tions, non-absorbed cyclo-synchrotron coupled 
to other radiation processes, acceleration pro- 
cesses (e.g. Fermi), pair creation/anihilation... 



4. First results 

Here we present a first use of the code. 
Several clues seem to indicate thermal plas- 
mas in AGN and in X-ray binaries although 
acceleration processes rather produce power 
law distributions. Simple particle-particle col- 
lisions by Coulomb interactions are too rare 
to account for this thermalization. Ghis ellini etl 
al. ( 1988] > suggested that exchange of energy 
between particles by exchanging synchrotron 
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